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HYPOCOERCIVITY FOR A LINEARIZED MULTI-SPECIES 

BOLTZMANN SYSTEM 

ESTHER S. DAUS, ANSGAR JUNGEL, GLEMENT MOUHOT, AND NICOLA ZAMPONI 


Abstract. A new coercivity estimate on the spectral gap of the linearized Boltzmann 
collision operator for multiple species is proved. The assumptions on the collision kernels 
include hard and Maxwellian potentials under Grad’s angular cut-off condition. Two 
proofs are given: a non-constructive one, based on the decomposition of the collision 
operator into a compact and a coercive part, and a constructive one, which exploits the 
“cross-effects” coming from collisions between different species and which yields explicit 
constants. Furthermore, the essential spectra of the linearized collision operator and 
the linearized Boltzmann operator are calculated. Based on the spectral-gap estimate, 
the exponential convergence towards global equilibrium with explicit rate is shown for 
solutions to the linearized multi-species Boltzmann system on the torus. The convergence 
is achieved by the interplay between the dissipative collision operator and the conservative 
transport operator and is proved by using the hypocoercivity method of Mouhot and 
Neumann. 


1. Introduction 

This paper is concerned with the proof of (explicit) spectral-gap estimates of the lin¬ 
earized Boltzmann operator for gas mixtures in the case of hard and Maxwellian potentials 
as well as the exponential decay of solutions to a multi-species Boltzmann system. Spectral- 
gap estimates and the large-time behavior of the mono-species Boltzmann equations were 
intensively studied in the literature, but are unknown for multi-species systems. First, we 
review the literature for the mono-species case. 

The study of the linearized collision operator, in the spatially homogeneous and hard- 
potential case, goes back to Hilbert [ 23 ] • For this operator, Carleman [5] proved the 
existence of a spectral gap. The results were extended by Grad HU for hard potentials 
with cut-off. Baranger and Mouhot [2| derived constructive estimates in the hard-sphere 
case. For Maxwell molecules, Fourier transform methods were employed in [37] to achieve 
explicit spectral properties. A spectral-gap estimate for the linearized Boltzmann operator, 
consisting of the sum of the linearized collision operator and the transport operator, was 
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first shown by Ukai [M]. Improved estimates (in smaller spaces of Sobolev type), still for 
hard potentials, were established in [ 22 ] • In [32], spectral-gap estimates for moderately 
soft potentials (withont angnlar cnt-off) were proved, improving and extending previons 
resnlts by Pao [SI]. Hypoelliptic estimates for the linearized operator withont cnt-off can 
be fonnd in [ 1 ] and references therein. A spectral analysis with relaxed tail decay and 
regnlarity conditions on the solntions was performed recently in an abstract framework 
[T9] . Dolbeault et ah [13] derived exponential decay rates in weighted spaces, which 
improves previons Sobolev estimates. For farther references, we refer to [301 Section 1.5]. 

Spectral properties of the linearized Boltzmann operator were already investigated by 
Grad [IB]. Based on these results, Schechter [32] located the essential spectrum of the 
classical collision operator in L^. The spectrum of the Boltzmann operator for hard spheres 
was also analyzed in U’ for p 7 ^ 2; see [23]. We refer to the recent work [13] for further 
results in W for 1 < p < cxd and more references. A detailed analysis of the resolvent 
and spectrum of the linearized Boltzmann operator can be found in [35l Section 2.2]. A 
complete analysis for the essential and discrete spectra for the linearized collision operator 
with hard potentials was performed in [28] . 

All these results are valid for the linearized mono-species collision operator. Our aim is 
to extend the spectral-gap analysis to the case of the linearized multi-species Boltzmann 
system modeling an ideal gas mixture. This is achieved by generalizing the coercivity 
method of [ 22 ], including quantitative estimates on the spectral gap for the multi-species 
collision operator. A crucial step of our analysis is the observation that the multi-species 
version of the iT-theorem implies conservation of mass for each species but conservation of 
momentum and energy only for the sum of all species. As a consequence, we need to study 
carefully the “cross-effects” of the collisions, i.e., how collisions between different species act 
on distribution functions which are elements of the nullspace of the mono-species collision 
operator. The crucial step is to relate these “cross-effects” to the differences of momentum 
and energy. Before stating the main results, we introduce the kinetic setting. 


1.1. The Boltzmann equation. The evolution of a dilute ideal gas composed of n > 2 
different species of chemically non-interacting mono-atomic particles (see m for chemically 
reacting gases) with the same particle mass can be modeled by the following system of 
Boltzmann equations, stated on the three-dimensional torus T^, 

(1) dtFi-^v-V^Fi = Qi{F), t>0, Fi{x,v,0) = Fi^i{x,v), (x, n) e x 

where 1 < i < n. The vector F = {Fi ,..., Fn) is the distribution function of the system, 
with Fi describing the ith species. The variables are the position x G T^, the velocity 
n G and the time f > 0. The right-hand side of the kinetic equation in ([1]) is the ith 
component of the nonlinear collision operator, dehned by 

n 

Qi{F) = '^Qij{Fi,Fj), l<i <n, 
i=i 
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where Qij models interactions between particles of the same {i = j) or of different species 
F^){v) = f Bij{\v - v%cos^){F^F'* - F,F*)dv*da, 

JR3x§2 

with the abbreviations F/ = Fi{v'), F* = Fi{v*), F'* = Fi{v'*), the three-dimensional nnit 
sphere S^, and 

, , , v + v* \v — v*\ v + v* In — 

(2) V = -h - -a, V = -- -a 

^ 2 2 2 2 

are the pre-collisional velocities depending on the post-collisional velocities {v,v*). These 
expressions follow from the fact that we assnme the collisions to be elastic, i.e., the mo- 
mentnm and kinetic energy are conserved on the microscopic level; 

11 11 

3 v+v =v + v, -n H —\v \ =-\v\ +-\v \. 

K J , 2 ' ' 2 ' ' 2 ' ' 2 ' ' 


The collision kernels Bij are nonnegative fnnctions of the modulus |n — n*| and the cosine 
of the deviation angle -d G [0, vr], dehned by cos'd = a ■ {v — v*)/\v — v*\. 

Although we will analyze a linearized version of Qi, let us recall the main properties of 
the nonlinear operator Q^. Using the techniques from m pp. 36-42], it is not difficult to 
see that Q := {Qi,..., Qn) conserves the mass of each species but only total momentum 
and energy, i.e. 



if and only if ip{v) G span{e‘l^\ ..., Uil, ^ 2 !, ^> 31 , |npl}, where is the Ah unit 
vector in R"' and 1 = (1,..., 1) G R”. It is shown in [11] that Q satishes a multi¬ 
species version of the F-theorem which implies that any local equilibrium, i.e. any function 
F being the maximum of the Boltzmann entropy, has the form of a local Maxwellian 

Mioc = (Mioc,i,..., Mioc,n) with 


Fi{x,v,t) = Mioc,i{x,v,t) 


(27r6>ioc(x,f))3/2 y 26>ioc(x,f) /’ 


where, introducing the total local density pioc = Pioc.i, 


Ploc^i 



^loc 


Plo 


n « 

^ / Fivdv, 


d\oc 



Fi\v — u\‘^dv 


are the (local) masses of the species, the total momentum and total energy, respectively. 

On the other hand, the global equilibrium, which is the unique stationary solution F to 
O, is given by M = (Mi,..., M„) with 


Fi{x,v) = Mi{v) 


Poo,i f 1^ ^ 00 1 \ 

(27rdoo)3/2 V 2 ^) 
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where now, setting = YJi=i Poo,u 

1 


Poo,i 


Fidxdv, Uoo = 


'T3xR3 


Poo Jt3xR3 


Fivdxdv, 600 = 


3po 


f I 

^i\v - u\ 

/T3xR3 



do not depend on (x, t). By translating and scaling the coordinate system, we may assume 
that Uoo = 0 and 6*00 = 1 such that the global equilibrium becomes 


(4) 


M,(n) 


Poo,i -|^|2/2 
(27r)3/2 


1 < i < n. 


1.2. Linearized Boltzmann collision operator. We assume that the distribution func- 
tion Fi is close to the global equilibrium such that we can write Fi = Mi + ' fi for some 

small perturbation fi, where Mj is given by (j4]). Then, dropping the small nonlinear 
remaining term, fi satishes the linearized equation 

(5) dtfi +V-Vxfi = Li{f), t>0, fi{x,v,0) = fi^i{x,v), (x, n) e T'^ x 

for I < i < n, where / = (/i ,..., fn) and the ith component of the linearized collision 
operator L = (Li,..., L„) is given by 

n 

Liif) = ^L,j{fi,fj), l<i<n, 

with 

fj) = M-'/' /,)) 

(6) = [ BijMl^^M*{h[ + h'*-hi-h*)dv*da, hr.= U 

2R3x§2 

Here, we have used M'*M'^ = M*Mj for any i, j, which follows from (|3]). Notice that we 
have chosen the linearization considered in, e.g., [29l|35]. Another linearization is given by 
Fi = Mi + MiPi (see, e.g., [27]), namely Lij{gi,gj) = M~^{Qij{Mi, Mipf) + Qij{Migi, Mf)). 
This choice gives the same results as with the linearization (I6|) since both linearizations 
correspond to the same space of solutions, but it turned out that the computations are 
easier using (]6|). 

The linearized Boltzmann system satishes an TT-theorem with the linearized entropy 

HU) = i fMv, 

= f‘UU)dv =: -(/,L(/))ti > 0, 

where (•, ■)i 2 is the scalar product on := We will prove in Lemma [5] that 

(/, L(/))i 2 = 0 if and only if fi lies in spanje^^^ ..., Xil, n 2 l, Xsl, |ppl}, which 

is the null space ff{L) of the linear operator L. The main aim of this paper is to show 
that, under suitable assumptions on the collision kernels, there exists a constant A > 0, 
which can be computed explicitly, such that for all suitable functions /, —(/, L(/))x 2 > 
where is the projection onto Mi^L) and "H is a subset of (see Theorem 
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[3] for the precise statement). This spectral-gap estimate, together with hypocoercivity 
techniques, allows us to conclude that exponential decay of the solutions f{t) towards the 
global equilibrium holds (see Theorem |lj) . 

1.3. Assumptions on the collision kernels. We impose the following assumptions on 
the collision kernels Bij arising in ([ 6 ]). 

(Al) The collision kernels satisfy 

Bij{\v — v*\, cos'd) = Bji{\v — v*\, cos'd) for 1 < i,j < n. 

(A2) The collision kernels decompose in the kinetic part > 0 and the angular part 
bij > 0 according to 

Bij{\v — n*|, cos'd) = — n*|) 6 ij(cos'd), I < i, j < 'n. 

(A3) For the kinetic part, there exist constants Ci, C 2 > 0, 7 G [0,1], and 6 G (0,1) such 
that for all 1 < i, j < 'u and r > 0 , 

Cir^ < <hij(r) <C 2 {r + r~^). 

(A4) For the angular part, there exist constants C 3 , C 4 > 0 such that for all 1 < f, j < 'u 
and'd G [0, tt], 

0 < bij{cos'd) < C 3 I slu'd! I cos'd!, %(cos'd) < C 4 . 

Furthermore, 

C’’ := min inf / min {bii{ai ■ as), bii(a 2 ■ a 3 )}da 3 > 0 . 

l<^<n (ji,cr2ES^ J^2 

(A5) For all 1 < 7 j < 'u, bij is even in [—1,1] and the mapping v 1 —)■ $h(!n!) on is 
locally integrable on and bounded as !'y! 00 . 

(A 6 ) There exists (3 > 0 such that for all I < i, j < n, s > 0, and a G [—1,1], we have 
Bij{s,a) < (3Bii{s,a). 

Following [22], since the functions bij are integrable, we dehne 
(7) := min / 6 jj(cosd) sinddd > 0. 

Jq 

Let us discuss these assumptions. The hrst hypothesis (Al) means that the collisions 
are micro-reversible. Assumption (A 2 ) is satished, for instance, for collision kernels de¬ 
rived from interaction potentials behaving like inverse-power laws. The lower bound in 
hypothesis (A3) includes power-law functions ^ij{r) = with 7 > 0 (hard potential) and 
7 = 0 (Maxwellian molecules). The assumption 7 > 0 is crucial since the linearized colli¬ 
sion operator in the mono-species case for soft potentials (7 < 0 ) with angular cut-off has 
no spectral gap [2]; however, degenerate spectral-gap estimates are possible [TBl [25] . The 
upper bound in (A3) means that the kinetic part is of restricted growth for both small and 
large values of !n — 'i;*!. In hypothesis (A4), the upper bound for bij implies Grad’s cut-off 
assumption. The positivity of in Assumption (A4) is used in the constructive proof 
of the multi-species spectral-gap estimate (Theorem [3]) via the 'mono-species spectral-gap 
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estimate which depends on see also the proofs of Theorem 1.1 in [2] and Theorem 6.1 
in [2B]. The positivity of is satisfied for the main physical case of a collision kernel 
satisfying Grad’s cnt-off, i.e. for hard spheres with — w*!, cos'd) = |n — n*!. Condi¬ 

tions (A1)-(A4) are also imposed in [21 [221127] for the linearized mono-species Boltzmann 
operator. 

Assnmption (A5) imposes technical conditions needed to verify the abstract hypotheses 
in [29]. More precisely, the evenness of bij is employed to show hypothesis (H2) (see section 
|5]) and the properties on are nsed to verify fl5T]) in hypothesis (HI). The conditions on 
are satished for hard and Maxwellian power-law potentials *hp(r) = with exponent 
7 G [0,1], for instance. Finally, condition (A 6 ) states that the ratio of the off-diagonal 
and diagonal collision kernels can be bonnded nniformly from above by a constant /3 > 0 . 
This hypothesis will be needed for the explicit compntation of the constants in Theorems 
[3] and 01 More precisely, (A 6 ) allows ns to estimate the mono-species part of the collision 
operator nsing the compntation of [26] ; see Lemma [TTl 


1.4. Notation and definitions. We call Dom(F’) the domain of an operator F and Im(/) 
the image of a mapping /. We introdnce the spaces Ll = L^(M^;M"'), = L^(T^ x 

M^;M"'), = H^{T^ x M^;M"'), and 


( 8 ) 


n=ue Li-. 


V=^feLi-. 

Here, Ui is the collision freqnency, given by 


n « 

< 

7^1 

n « 

• 1 Jrs 


CXO 


(9) 


n « 

/ 

j=l 


’XS2 


— n|, cos7})M*dv*da, i = 1 ,... ,n, 


For Maxwellian modelcules 4’p(r) = const., the collision freqnency is constant bnt for 
strictly hard potentials ‘hp(r) = with 0 < 7 < 1, z/j is nnbonnded. In fact, it satishes 
z/o(l -|- |n|)'’' < i'i{v) < z/i(l -|- Inl)"^ for some constants > zzq > 0 [29l p. 991]. In the 
physically most relevant case of hard spheres (7 = 1 , % = 1 ), the collision freqnency can 
be compnted explicitly, see formnla (2.13) in [TUI Section 7.2]. For more properties of the 
collision freqnencies, we refer to 0 Section HI.3]. If the collision freqnencies are bonnded, 
LL = Li,. Generally, z/j is nnbonnded and so, "H is a proper subset of Li. The norm on Li 
(and similarly for the other spaces) is defined by 


ll/ll 


2 

Ll 



for / = (/i,...,/„) G Li. 


We distinguish the following linear operators. We dehne the operator A = (Ai,..., A„) : 
Dom(A) —)■ Li by 


Mf) = Ffi, i = l,...,n, 
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where Dom(A) = {/ G : A/ G L^} = V. It is closed, densely defined, selfadjoint and, by 
Lemma [7] below, coercive. The linearized collision operator L : Dom(L) ^ L^, introduced 
in section [L2l can be written as L = K — A, where K := L + A, or, more explicitly, 

n „ 

Ki{f) = V / + h'* - h*)dv*da, i = 1 ,..., n. 

jR3xS2 

It was shown in [1] that K is a compact operator in L^. Thus, Dom(L) = Dom(A) = V 
and L is closed and densely defined. Furthermore, L is nonpositive and selfadjoint on L^. 
We define the transport operator 

T = v - Vx ■ Dom(T) Ll, 

where Dom(T) = {/ G ; v-A/xf G Finally, we consider the linearized Boltzmann 

operator 

B = L-T: Dom{B) Ll, 

which is unbounded, closed, and densely defined with Dom(i?) = Dom(L) fl Dom(T). 

We denote by Afi^A) and 'R.{A) the kernel and range of a linear operator A, respectively. 
Its resolvent set is denoted by p{A) and its spectrum by cr(A) = C\p(A). For a linear 
unbounded operator A with a{A) C (—oo,0], we say that A has a spectral gap when the 
distance between 0 and (t(A)\{ 0} is positive. Finally, the essential spectrum of A is defined 
as the set of all complex numbers A G C such that A — AJ is not Fredholm, where I is the 
identity operator. We refer to section [3] for details regarding this definition. 

1.5. Main results. In this subsection, we state the main results of the paper. The first 
result is a geometric property of the essential spectrum of the linearized collision operator 
L and the linearized Boltzmann operator B = L — T. 

Theorem 1 (Essential spectrum of L and L — T). Let the collision kernels Bij satisfy as¬ 
sumptions (Al)-(Af) and set J = UjL;^Im(r'i) C [z/o,cxo), where oq = minj=i^...^„ sup^g^s i^i{v) 
> 0 (see Lemma Then 

(Tess(T) = -J, a^L - T) = {A G C : 3?(A) G - J}. 

Remark 2. We observe that if hm|^|_,.oo r>i{v) = oo for i = 1,... ,n then 

(Tess(T) = (-00, -Z/q], (Tess{L - T) = {A G C : 3f?(A) < -Uq}. 

Indeed, under the assumption z/j(n) —)■ oo as |n| —)■ oo, the continuity of z/j, and the 
Weierstrafi theorem show that J = [z/q, oo). Thus, the essential spectrum of the linearized 
multi-species collision operator is very similar to the mono-species operator, where z/q cor¬ 
responds to the infimum in of the single collision frequency; see m Section 3] and [2H1 
Prop. 3.1]. □ 

The proof of Theorem [1] is based on perturbation theory [231 Chap. IV] and is similar to 
the proof for the mono-species collision operator [35], but we show new explicit spectral-gap 
estimates related to the particular structure of the kernel in the multi-species case. More 
precisely, we write L = K — A as described in section 11.41 It turns out that K = L + A 
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is compact on Ll (see section [3] for details). Weyl’s theorem [20l Theorem S] states that 
the essential spectrum oi L = K — A coincides with that of —A. Thus it remains to show 
that dess (A) = J. This is done by using Weyl’s singular sequences, which allow for a 
sufficient and necessary condition for A G C being an element of the essential spectrum of 
the selfadjoint operator A. 

The proof of the second statement in Theorem [1] is more involved since K is not compact 
on Ll ^ and hence, Weyl’s theorem cannot be applied directly. The idea is to employ an 
extended Weyl theorem, which states that the essential spectrum is conserved under a 
relatively compact perturbation pTl Section IV.5.6, Theorem 5.35]. Indeed, if K is relatively 
compact with respect to A + T then (Jess(A — T) = aess{K — (A + T)) = —cregg(A + T), and 
it remains to compute the essential spectrum of A + T. 

The next theorem concerns an explicit spectral-gap estimate. It is the main result of the 
paper. 


Theorem 3 (Explicit spectral-gap estimate). Let the collision kernels B^j satisfy assump¬ 
tions (Al)-(Af). Then there exists a constant A > 0 such that 

(10) -(/,i(/))L;>A||/-ll‘(/)||?, forallfev, 

where 11 ^ is the projection onto the null space M{L). If additionally hypothesis (A6) holds, 
the constant A can be computed explicitly: 


A 


rjD^ 

SCk 


7] = mm 


’ IQCk + 


where C™, , and Ck are defined in fl38l) . fH3l) . and fll5l) . respectively. 


Note that the constant C™ depends on the mono-species spectral-gap constant via 
fl38ll below. We present two proofs of this theorem. The hrst proof is non-constructive 
and relies on an abstract functional theoretical argument, based on the decomposition 
L = K — A and Weyl’s perturbation theorem. This abstract spectral-gap estimate is 
proved in Lemma dU) The second proof provides a constructive spectral-gap estimate, 
generalizing the result in [27] (also see [251 Theorem 6.1]) from the mono-species to the 
multi-species case. For this, we split the operator L = L™ -|- L^ in the mono-species part 
L™ = (L™,..., L™) and the bi-species part Lf = {L \,..., L^), 


rr(/i) = /i), r?(/) = V *)■ 


The proof consists of four main steps. 

Step 1: Coercivity of the mono-species operator L^. The bi-species part of L satishes 
— (/, L^(/))i 2 > 0 for all f eV. Furthermore, the results of [25] Theorem 6.1] show that 
for the mono-species part, 

(11) - (/, L'-U))lI > C^Wf - n^U)\\li for / G Dom(L”^), 

where the constant > 0 can be computed explicitly and il™ is the projection onto 
M{L'^) (see Lemma ITT]) . Inequality flTT]) may be interpreted as a coercivity estimate for 
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L™ on It is related to the “microscopic coercivity” in [131 Section 1.3] for the 

mono-species setting. Hence, we obtain the “naive” spectral-gap estimate 

-(/,^(/))Lg > C^Wf - n^ifWn for fev. 

This estimate is not sharp enongh for the mnlti-species case since we need an ineqnality for 
all / — n^{f) G A/'(T)-‘- and not only for / — n'^{f) G C By projecting 

onto only, we neglect the “cross-effects” coming from the bi-species part of the 

collision operator. Thns, we need a better estimate for —(/, L(/))^2, which is achieved as 
follows. 

Step 2: Absorption of the orthogonal parts. The contribntion f'^ '■= f — n^{f) in 
— (/, Lf{f))L 2 can be absorbed by the "H norm of /■'■ (see Lemma [12]), giving for a certain 
h > 0, 

-(/. mu £ (c” - mnwii - Imu), 

Step 3: Coercivity of the hi-species operator Lf. The projection n"^{f) depends on the 
velocities Ui and energies e* of the ith species, and thns, the cross terms can be bonnded 
by the differences of momentnm and differences of energies, 

n 

*j=i 

for some constant C > 0. This is the key step of the proof. The inequality may be con¬ 
sidered as a coercivity estimate for the bi-species operator. A key observation is that the 
differences of momenta and energies converge to zero as / approaches the global equilib¬ 
rium. Whereas fITTD acts on flT^ gives an estimate on the orthogonal complement 

Step f: Lower bound for the differences of momenta and energy. The last step consists in 
estimating the differences In,- — I and le,- — I from below by the error made by proiecting 
onto instead of AA(L)^: 

n 

(If, - up + le, - ep) >C\\f- nmfn - 2C-||/ - 17'“(/)||L 

Putting together the above inequalities. Theorem [3] follows; we refer to section 0] for details. 

As a consequence of the spectral-gap estimate, we are able to prove the exponential 
decay of the solution f{t) to ([5]) to the global equilibrium with an explicit decay rate. 

Theorem 4 (Convergence to equilibrium). Let the collision kernels Bij satisfy assumptions 
(A1)-(A5) and let fj G Hf.^. Then the linearized Boltzmann operator B = L — T generates 
a strongly continuous semigroup on which satisfies 

(13) ||e^^(J-77'')bi„ <Ce-”*, f > 0, 

for some constants C, r > 0. In particular, the solution f{t) = fi to ([5]) satisfies 
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where foo ■= is the global equilibrium of ([5]). Moreover, under the additional 

assumption (A6) and lower bound in (Af), the constants C and r depend only on the 
constants appearing in hypotheses (H1)-(H3) in section\E and in particular on A defined in 
Theorem O 

The idea of the proof is to employ the hypocoercivity of the linearized Boltzmann op¬ 
erator L — T, using the interplay between the degenerate-dissipative properties of L and 
the conservative properties of T. The aim is to hnd a functional G[f] which is equivalent 
to the square of the norm of a Banach space (here, 



for / G 


leading to 
(15) 



where ai, K 2 , n > 0 and f{t) = fj. These two estimates yield exponential convergence 
of f(t) in Hl^. It turns out that the obvious choice G[f] = Ci||/||i2^ + C 2 \\'^xf\\Li„ + 
CsIIV^/IIlz^ does not lead to a closed estimate. The key idea, inspired from [SB] and 
worked out in is to add the “mixed term” C 4 (Vxf,'^vf)Ll„ to the dehnition of G[f]. 
Then 



and the last term can be estimated in terms of expressions arising from the time derivative 
of the other norms in G[f(t)]. Thus, choosing q > 0 in a suitable way, one may conclude 
that ffT5|) holds. 

In [29], the calculation of flT5|) is reduced to the validity of certain abstract conditions 
on the operators K and A (see section [5]). These conditions state that A is coercive in a 
certain sense, K has a regularizing effect, and L = K — A has a local spectral gap. The last 
condition is proved in Theorem jS] while the other conditions follow from direct calculations, 
since the operators K and A are given explicitly. As a consequence, the proof of Theorem 
m essentially consists in verifying the abstract conditions stated in [29]. In contrast to the 
estimate of Theorem [31 where the multi-species character plays a role in the spectral-gap 
estimate, there are no “cross-effects” here and the same modihed functional G[f] as above, 
including the mixed term, can be used. However, the decay rate r changes, since the 
constant in hypothesis (H3) (see section [5|) differs in the mono- and multi-species case and 
r depends also on that constant. 

We hnish the introduction by commenting possible generalizations. First, the conver¬ 
gence result based on hypocoercivity requires some regularity on the initial data, namely 
fi G The extension of the exponential decay to initial data from ^ might be done 
by using the method of (TB), which is based on a high-order factorization argument on the 
resolvents and semigroups. The proof of exponential decay is expected to be constructive 
and to preserve the optimal rate. 

Second, it seems to be not trivial to extend the results to the whole-space case. The 
problem is that one loses the compactness in the x-space. One possibility is to assume some 
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confinement potential which, nnder some appropriate weighted Poincare inequality, can 
yield compactness of the resolvent and hence a spectral gap. For instance, Duan [T3] used 
non-constructive techniques to prove decay rates for the mono-species linearized Boltzmann 
equation. Still in the mono-species case, with one-dimensional collisional invariants and 
using constructive methods, the decay is investigated by, e.g., Herau and Nier [22] and 
Villani [36], working in the space Hi and by Herau [21] and Dolbeault, Mouhot, Schmeiser 
[T3] . working in the space The tasks in the multi-species case are first to extend the 
non-constructive methods, which probably does not contain new difficulties, and second to 
devise a constructive method, which is more involved and work in progress; see [T2] . 

Third, a Cauchy theory for the full nonlinear multi-species Boltzmann equation in a 
perturbative regime is work in progress [7]. 

The paper is organized as follows. In section |2l some properties of the linearized collision 
operator ([6]) are collected. Theorem [1] on the essential spectrum of L and L — T and the 
abstract spectral-gap estimate in Theorem [3] are proved in section (3] We present a second 
proof of Theorem [3] in section 01 by exploiting the conservation properties and leading to 
explicit constants. Finally, Theorem 0] is shown in section 01 

2. Properties of the kinetic model 

We show some properties of the linearized collision operator ([6]) and the collision fre¬ 
quencies 0. Let assumptions (A1)-(A4) hold. First we prove an FT-theorem for ([6]). 

Lemma 5 (FT-theorem for the linearized collision operator). It holds {f,L{f))i 2 < 0 for 
all f eV and (/, L{f))i 2 = 0 if and only if f E J\r{L), where 

Af{L) = [f E LI : 3q!i, ..., an, e G M, u E VI < i < n, 

fi = Ml^‘^{ai + U-V + e\v\‘^)}, 

and Mi is given by (0]) . 

The proof is similar to the mono-species case except that the elements of the null space 
of F depend on the total mean velocity u and total energy e instead of the individual 
velocities and energies. Therefore, we give a complete proof. We note that an Ff-theorem 
for the nonlinear Boltzmann operator for a mixture of reactive gases was proved in [m. 

Proof. By the change of variables {v, v*) ha {v*, v) and {v, v*) ha (n', v'*) and the symmetry 
of Bij (assumption (Al)), we can write for f E L^, 

1 n „ 

(/, Lif))L 2 = -T V / B,,MiM;{h[ + hi* - F, - h*fdv*dvda, 

where we recall that hi = M- This shows that {f,L{f))L 2 < 0 for all f E V. 

Moreover, (/,F(/))i 2 = 0 if and only if 

(16) h'i + h'* — hi — h* = for all (n, v*) G x M^, 1 < i, j < n. 
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It is shown in m pp. 36-42] that flT6|) for i = j implies that hi has the form hi{y) = 
ai + Ui ■ V + ei\v\^ for snitable constants aj, Cj G M and Ui G Inserting this expression 
into (fT6l) leads to 

(17) Ui ■ {v' - n) Uj ■ {v'* - V*) + ej(|n'|^ - |n|^) ej{\v'*\‘^ - |n*|^) = 0 

for 1 < i,j < n. We consider the particular type of collisions with v' = v*, v'* = v, and 
|n| = |n'|. For such collisions, a = (v* — v)/\v* — v\. Then the above equation becomes 

{ui — Uj) ■ {v' — v) = 0, 1 < i,j < n. 

By rotating the velocities n, n' in all possible ways, we deduce that {ui — Uj) ■ w = 0 for all 
tc G and thus, Ui = uj for all 1 < i, j < n. We set u := ui. This fact, together with the 
conservation of momentum v' — v + v'* — v* = 0, implies that flTTI) becomes 

ei(|n'p — |n|^) + ej{\v'*\‘^ — |n*|^) = 0, 1 < f, j < n. 

Taking into account the conservation of energy = 0, we infer that 

(cj — ej)(|n'p — |np) = 0 and consequently, e, = ej for all 1 < i,j < n. Set e := Ci. We 
have shown that (/, L{f))i ,2 = 0 if and only if there exist ai,..., e G M and m G 

such that for 1 < ^ < n, fi{v) = Ml^'^[ai + u- v + e|np). These functions clearly belong to 
AA(L), which hnishes the proof. □ 


The next result is concerned with the stationary solutions of (jS]). 


Lemma 6. The global equilibrium foo = (/oo,i, • • •, foo,n) of ([5]), i.e. the unique stationary 
solution, is given by 

/oo,i(y) = Ml^‘^{ai + U-V + e|n|^), I < i < n, 

where e G M and n G are uniquely determined by the global conservation laws of 
mass, momentum, and energy, i.e. by the equations 


[ fij)'ip{v)dv 

Jr3 


0 , I < i < n, 


for fj{v) = l,Vi,V 2 ,V 3 , |np, where fjj are the initial data. 

Proof. First, we claim that M{B) = M{L) n M{T), where B = L — T and T = v ■ Va, are 
considered on x The inclusion Af{L) r\Af(T) C J^{B) being trivial, let / G J\f{B). 
Then, using the skew-symmetry of T, 


0 = (/,B(/))i|.. 






Lemma E] shows that / G N'{L). But this implies that T(/) = L{f) — B{f) = 0 and hence 
/ G Af(T). This shows the claim. Let foo be a stationary solution. Then foo G J^{B) and 
by our claim, / G M^L) PMiT). Since J\f{T) = {/ G : Vxf = 0} [HI Lemma B.2], foo 
does not depend on x. Because of foo G Lemma [5] shows the result. □ 


Finally, we prove that the collision frequencies ([9]) are strictly positive with bounded 
derivative. 
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Lemma 7. Let Assumptions (A2)-(A4) hold. The collision frequencies (E]) satisfy 

(18) min inf Ui{v) > := 23 ^/ 2 > 0 , 

y'TT \ 2 J 

where Ci > is given by assumption (AS), > 0 is defined in ([7]), poo := YTj=iPj,oo 

(see ^), and V is the Gamma function. Furthermore, if additionally (A5) holds, then 
G implying that \ni{v)\ < C,y{l + |w|) for some C,y > 0 and for all v G R^, 

i = 1,... ,n. 


Proof. The decomposition of Bij, according to assumption (A2), implies that 


n « 

/ ^ij{\v — v*\)M*dv* I bij{cosd)da. 

j=i 


/S 2 


The integral 


Cij := / bij (cos ■&)da = 271 / by {cos d) sin Add 

is 2 Jo 

does not depend on v or v*. We conclude from (A2)-(A4) that 
(19) 


z/i(u) = (27r) ^^^^CijPooj 
i=i 

Cl^Voo 


> 


(27r)3/2 




Observe that the function 


G(u) := / |u-u*|T'e-h*lV2^y* 


is uniformly positive since the transformation v* i—)■ —v* and the elementary inequality 

|u — U*|3' + |u + V*\"' > |(u — V*) + (u + u*)|'^ = 2'^|u*|'’' 

for 7 G [0,1] lead to 
1 


= o / ~ + |u + u*|'^)e ^‘^dv* > 2'^ ^ 

2 Jm3 


|y*|7e-KI Gdv* = 2^-^G{0). 
Actually, using spherical coordinates and the change of unknowns s = 'r^/2, 

G( 0 ) = 471 • 


Inserting the above estimate on G{v) into flT^ shows flTSl) . 

It remains to prove that G L“(K^)- To this end, we compute 


|Vp,(u)| = (27r)-3/2 


^ ^ GjPi 


oo,j 


4>d(|77 -U*|) • 


V — V 


\v — V 
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n „ 

Jv 

For given i? > 0, we decompose 

l^pd^DI = l^p(kl)IX{H<ii}(w) + 

Assumption (A5) means that there exists i? > 0 such that 

l-i-td • DIxii kr) 6 i;(Rt and |4p| ■ |)|xih>h) 6 ir(R")- 

Thus, the right-hand side of (l20ll is bounded since it can be written as the sum of two 
terms, each of which is the convolution of an and an L°° function. This shows that 

e □ 

Remark 8 . We observe that z/j is generally not bounded since the kinetic part *hp(r) may 
grow like r as r —)■ oo. It is possible to show that z/* is bounded if ^ij is bounded. The 
unboundedness of Pi implies that the spaces and "H are not isomorphic. □ 


3. Geometric properties of the spectrum 


In this section, we prove Theorem [U and the spectral-gap estimate flTOll in Theorem |3] by 
using arguments from functional analysis. 

First, we study the essential spectrum of L and L — T. There exist several definitions 
of the essential spectrum of a linear operator. Given a linear, closed and densely defined 
operator A : Dom(A) C X —)■ X on a Banach space X, we define 

c’‘ess(A) = {AgC:A — A/is not Fredholm}. 

We recall that a linear, closed, and densely defined operator A is Fredholm if its range TZ{A) 
is closed and both its kernel and cokernel are finite-dimensional. For other definitions of 
the essential spectrum, we refer to |20]. The essential spectrum is closed and conserved 
under compact perturbations, i.e., the bounded operators A and B have the same essential 
spectrum if A — R is compact (Weyl’s theorem; see [201 Theorem S]. 

If X is a Hilbert space and A is selfadjoint, it holds cress(A) C M and for given A G 
M, we have A G dess (A) if and only if A — AJ is not closed or the kernel of A — AJ is 
inhnite dimensional. (This follows from the fact that TZ{A — A/)-*- = J\f{A — XI) for 
closed, selfadjoint operators A [2ll Chap. V.3.1].) Moreover, Weyl’s criterion holds [3^ 
Lemma 5.17]: A G (Tegs(A) if and only if A — AJ admits a singular sequence, i.e. a sequence 
(/fc) C Dom(A) such that (i) ||/fc||x = 1 for all fc G N; (ii) ||(A — A/)/fc||x —)■ 0 as fc —)■ oo; 
and (hi) {fk) has no convergent subsequences in X. 

We decompose L as L = K — A, where K = (Xi,..., X„), A = (Ai,..., A„), and 



+ h'* 


Ai(/) = Uifi, I <i<n, 


h*)dv*da, 


(21) 
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and the collision frequencies Vi are defined in ([9]). We recall from Lemma [3 that they satisfy 
Pi('i') > pq > 0 and \vi{v)\ < + |n|) for alH = 1,...,n and v G 

Proof of TheoremUl Since K is compact on X = [H Prop. 2], it follows that aess{L) = 
(Jess(—A) = —cress(A). Thus, wc will hrst study the essential spectrum of A. The proof is 
divided into several steps. Recall that J = U”^]^Im(z/j) C [uq, oo). 

Step 1: J C (Tess(A). Let A G J. There exists j G n} and n G such that 

A = Pj (n). We dehne the sequence {fk) C "D by 

fk,i{v) = (27r(Tfc)"^/‘^exp ^ if * = 2) fkA'^) =0 if * 7^ J, 

where ak = l//c, k G N. Clearly, condition (i) for the singular sequence is satished. 
Furthermore, 

n p 

II (A - A/)/fc||i2 = / {^i{v) - XffkAvfdv 

i=i 

= ( 27 rcrfe)"^A f _ u.{v)Yexp ) dv. 

7r3 V / 

The limit of a sequence of Gaussians with variance tending to zero converges to the delta 
distribution (in the sense of distributions), which means that 

( 27 r(Tfc)“^A f exp f ———^ j dv —)> u{v) as A; —)■ oo 

Jr 3 \ 2ak J 

for all functions u G with polynomial growth at inhnity. Since Wi{v)\ < C^(l + |n|), 

this condition is satished and we conclude that ||(A — A/)/fc||L 2 —)■ 0 as /c —)■ oo, showing 
that condition (ii) holds. 

Let ns assume by contradiction that condition (iii) does not hold. Then there exists a 
subsequence (/^J of (/*,) that converges in to some function f ^ Li. As a consequence. 
I/p in L\ as I ^ oo. In particular, / G L^. However, the distributional limit 
dv and the uniqueness of the limit imply that 5^ = |/P G L),, which is absurd. 
Thus, condition (iii) holds, and we infer that A G cress (A). Then, since dess (A) is closed, 
J (Z cregg(A). 

Step 2: cregg(A) C J. Let A G M\J. Then there exists a constant c > 0 such that for all 
n G and i = 1,... ,n, |r'i(n) — A| > c. If {fk) C V with ||/fc||L2 = 1 for all A; G M, we 
have 

n p 

iivfdv > V] / fkA'^fdv = > 0 

for all A: G M. Thus, condition (ii) cannot hold which implies that A ^ dess (A). 

Steps 1 and 2 imply that dess (A) = J■ 

Step 3: {A G C : 3f?(A) G J} C dess (A + T). Let A G C be such that 3f?(A) G J. It 
follows from Step 1 that 3f?(A) G dess (A). Since A is selfadjoint on the Hilbert space Li, 
A — 3?(A)/ is not closed or the kernel of A — 3?(A)/ is inhnite dimensional. As the operator 


||(A-A/)M|i2,= 


n p 

- Xffk, 
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A — 3?(A)/ is closed, its kernel must be infinite dimensional. Therefore, there exists a 
sequence (/*,) C such that A(/fc) — 3fJ(A)/fc = 0 and (/fc, fi)Ll = for k, (. e'H. Let us 
define (j){x,v) = exp(i9(A)a; • n/|np) and = 4>fk £ Since |0| = 1, we have 

(22) {gk, gi)Li^^ = {fk, f^Li = ^ki for fc, ^ G N. 

Furthermore, (j) G Dom(T) and T(0) = i9'(A)0 for n 7 ^ 0, and thus, 

(A + T - XI)gk = 0(A - ^{X)I)fk + fk{T - i3(A)/)0 = 0, 

which shows that gk G J\f{A + T — XI) for A; G N. This fact, together with relation ([22]), 
implies that J\f{A + T — XI) is infinite dimensional. As a consequence, A + T — AJ is not 
Fredholm and A G dess (A + T), which proves the claim. 

Step 4- {A G C : 3?(A) ^ J} C p(A + T). Clearly, this gives 

(Tess(A + T) C a(A + T) C {A G C : 3?(A) G 7}. 

Let A G C be such that )R(A) G M\J. We show first that J\f{A + T — XI) = {0}. We 
assume by contradiction that there exists / G Dom(A + T) satisfying ||/||l 2 ^ > 0 and 
(A+T —A/)/ = 0. In particular, there is an index £ G {1,..., n} such that ffdxdv > 

0. Then, multiplying = Xfi by (the complex conjugate of fi) and integrating 

in X M^, we obtain 


(23) 



I'ilfil'^dxdv + 


K3 JT3 



f■ Vxfidxdv = X 


K3 Jt3 



Ifil'^dxdv. 


R3 JT3 


By the divergence theorem, the real part of the second integral vanishes. 


(24) 


23? 



■ Vxfedxdv 



{fiv ■ '^xfi + fiv ■ Vxfi)dxdv 
r’ ■ A/x\fe\'^dxdv = 0. 


Then, taking the real part of (|23|) . we infer that 

^ / r 3 Its Mfi?dxdv 
/ r 3 fj3 if^i^dxdv 


Consequently, inf^s < 3?(A) < supj^s and, thanks to the continuity of 3?(A) G 
Im(z/£) C J, which is a contradiction. Thus, Af(A + T — XI) = {0}. Similarly, we can show 
that M{{A + T- XI)*) = M{A + T*- XI) = {0} as well. 

The operator L — T is closed [SS] Theorems 2.2.1 and 2.2.2]. Thus, the boundedness 
of the compact operator K and the stability of closedness under bounded perturbations 
[241 Chap. Ill, Problem 5.6] imply that A + T = K — {L — T) is closed (and also densely 
defined). Hence, Il{A + T — XI)^ = A/'((A + T — XI)*) = {0}, meaning that A + T — AJ is 
invertible. If / G is given, there exists u G Dom(A + T) such that (A + T — XI)u = /, 
which translates into 


(25) {uj - 3?(A))mj + (T - i‘A{X))uj = /, j = 
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We point out that, since Vj is continuous, lm{uj) is an interval (or a point, in case that uj 
is constant). This fact and the assumption 3?(A) ^ J imply that either Uj — 3fJ(A) > 0 in 
or Uj — 3fJ(A) < 0 in This means that the sign Sj of uj — 3fJ(A) is constant in for 
j = 1,... ,n. By multiplying fl2S]) by SjUj, integrating over x taking the real part, 
and summing over j = 1 ,..., n, we find that 





Jm.^ Jt^ 


Wj-^{X)\\uj\‘^dxdv = 2 

i=i 



M3 JT3 


Ujfj + Ujfj)dxdv. 


The (real part of the) second term in fl25l) vanishes after integration; see fl2Tll . Since 
A G M\J, by definition of J, there exists ca > 0 such that \uj — 3fJ(A)| > cx in for all 
j = 1,... ,n. Then the Cauchy-Schwarz inequality shows that UmIIlj^ < Ca^||/||l2 This 
means that {A + T — is bounded, so A G p{A + T). 

Steps 3 and 4 show that (Jess(A -|- T) = {A G C : 3fJ(A) G J}. 

Step 5: cress(—A — T) = a^ss^K — A —T). The operator K is compact on L‘1 but not on 
so the claim does not follow from the original form of Weyl’s theorem. Instead we 
will employ the fact that the essential spectrum is conserved under a relatively compact 
perturbation [2U Section IV.5.6, Theorem 5.35]. More precisely, we prove that K is rela¬ 
tively compact with respect to A -|- T, i.e., := K{A + T — zl)~^ is compact on for 

some z G C with 3fJ(z) G M\J. (Notice that by Step 4, 2 ; G p(A -|- T).) Then 


(26) aess(i^ - A - T) = aess(-A - T) = {A G C : 3f?(A) G - J}. 


The second identity is a consequence of Steps 3 and 4. 

To prove the compactness of we introduce the space W := of sequences 

/ = {fm) C Ll with the canonical norm \\f\\w = ll/mllij)^/^- Clearly, W is 

a Hilbert space with the scalar product {f,g)w = 9m)Ll- Furthermore, we 

introduce the Fourier mapping F : L^^(T^ x M^) ^ PF by 


F{f) = {fm), fm{v)= [ e-^-^^-^f{x,v)dx for m G T G RI 

At3 

This mapping is bounded, invertible, and has a bounded inverse. We wish to show that 
B^ = FBzF~^ : IF —)■ IF is compact. Then also B^ = F~^BzF is compact as a composition 
of a compact and two bounded operators. This idea is due to Ukai; see e.g. [321 Section 
2 . 2 . 1 ]. 

Since K and A do not depend on x, it holds that B^ = K{A+T—z)~^, where T = 27rin-m. 
Let (/^^^) = {fm') C IF be a bounded sequence in IF, i.e., there exists cq > 0 such that 
for all /c G N, 


ll/“'ll?,- = E S "O' 

meZ3 


(27) 
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As G M\J, there is a constant > 0 snch that for all i = 1 ,... ,n and v G 

\i^i{v) + 27rin ■ m — z\ > c^. Thus, 


n « 

||(A + f - = 5^ / \{iz,{v) + 27rto . m - 

i=l 

n « 

<c-E/ l/Sr<i<' = c"ll/.'t'lli;- 

i.l ■'"* 


Summing these inequalities over m G we infer that 

||(A + f - < c«cf. 


Consequently, the sequence := {A + T — z) is bounded in W. In particular, for 

any s G Z^, || 5 f ?^||^2 < Emezs < <^oC~^- Hence, for any s G Z^, the 

sequence C is bounded in L^. Since K : ^ is compact and Z^ is countable, 

we may apply Cantor’s diagonal argument to hnd a subsequence of {g^^'^) such that 

{K{gm‘'^)) is convergent in as £ — )> oo, for all m G Z^. 

We will show that is a Cauchy sequence in W. To this end, let i, s, N ^ N. 

We write 


(28) ’) - K(st^)\\k 

mez® 

= E \\i<(a'y)-K(at’)\\h+ E IDlsS't-AlsS-tlll;. 

|m|<A |m|>A 


where \m\ = Ef=i ^ First, we consider the second sum on the right-hand 

side. Denote by || • \\^{lI) the norm in the space of linear bounded operators on L^. By 

(k) 

the dehnition of gm ■, we obtain 

(29) E IA(9S'’) - K(g>J:->)\\% = E + 2 ™ . m - 2)-‘(/t*'> - ft')\\li 

|m|>A |m|>A 

|m|>A 

For the operator norm, we employ Prop. 2.2.6 in [33], which can be applied since 3fJ(2;) G 

M\7: 

||A'(A -I- 27rin • m — < Ci(l -|- |m|)““ for all m G Z^ 

for some suitable constant ci > 0 (depending on z) and a suitable exponent a G (0,1) 
(actually, a = 4/13). Let 0 < (3 < 2a/3. By Holder’s inequality and fl27)) . we estimate 

\\K{A + 27riv-m-z)-^\\%(^Li)\\fm\\l 2 <Co^ci Y (1 + 

|m|>A^ |m|>A^ 
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1-/3/2 


Using this estimate in (1^ . it follows that 

sup < 2C0Cl 



/3/2 


The choice of 13 implies that 2a/(3 > 3 and hence, the sum over \m\ > iV is hnite. In 
particular, X]|m|>N(^ + —)■ 0 as —?■ cx). Asa consequence, for given e > 0, there 

exists G N such that 

sup E 

Finally, since {K{gia^ )) is convergent in for all m G there is a number rj = 77 (e) > 0 
such that for all £, s > r], 

E IIM9S'’)-M9S-’)lli;<|. 

\m\<Ne 

Thus, choosing A^ = Ag in fl28|) . we deduce that is a Cauchy sequence in the 

Hilbert space W and consequently, it is convergent. This shows that : IT ^ lU is a 
compact operator and (l26ll holds. This hnishes the proof of Theorem [H □ 


Next, we show the spectral-gap estimate for the linearized collision operator L = K — A, 
i.e. the hrst statement of Theorem |3l Since K is compact on L^, it remains to prove that 
A-.VcLl^Ll is coercive. 

Lemma 9. Let (Al)-(A4) hold. Then the embedding TL ^ is continuous and A : P — )■ 
Ll, defined in CD, is a linear unbounded operator with the property 

(30) {f.m)Li = \\f\?n>C\\f\\l. for fen 

for some C > 0. Moreover, A can be extended by density to a linear bounded operator 
A ; 'H —)■ H', where TL' is the dual ofTL with respect to the scalar product. In particular, 
the mapping "H —)■ M, / 1 —)> (A(/), /) is continuous, where (•, •) denotes the duality pairing 
between TL' and TL. 


Proof. The strict positivity of z/j in (see Lemma [7]) implies that the embedding TL ^ 
is continuous. Then the dehnitions of Aj and TL show that for all f ^ TL, fl30|) holds. For 
given f eTL, the element A(/) = (/i^i,..., fn^n) can be identihed with the linear bounded 
operator TL M., g ^ X]r=i Irs difi'^id'v and consequently, A(/) G TL'. It is immediate to 
see that ||A(/)||^/ = H/H-H; so that A : "H —)■ "H' is isometric and thus bounded. Moreover, 
it follows that "H —)■ M, / e-)■ (A(/), /), is continuous. □ 
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The following result provides a spectral gap for general operators which decompose into 
a compact and a coercive part. 

Lemma 10. Let Hq and LL be Hilbert spaces such that H ^ Hq continuously and let 
L ; n H' be a linear bounded operator such that L = K — A with linear bounded 
operators A ■. H ^ H,' and K : Hq Hq. Furthermore, assume that 

(i) for all f eH, {L{f), /) < 0 with equality holding if and only if f & Af{L); 

(ii) the operator K : Ho —?■ Ho is compact; 

(iii) there exists Co > 0 such that for all f eH, (A(/), /) > Co||/|||^. 

Then there exists a constant Ci > 0 such that 

-W), /) > cml for anf€H n V(L)^. 

Proof. We argue by contradiction. Let (/„) C HnAf{L)-^ be a sequence such that H/nll-H = 
1 for n > 1 but {L{fn), fn) ^ 0 as n ^ cxo. Since (/„) is bounded in the Hilbert space H, 
there exists a subsequence, which is not relabeled, such that fn^f weakly in H. Because 
of the continuous embedding H Ho, also fn^f weakly in Ho- Since fn G Af{L)-^ and 
Af{L)-^ is weakly closed by Mazur’s lemma, / G Af{L)-^. As the operator K : Ho —)■ Ho is 
compact, by hypothesis (ii), the weak convergence of (/„) in Ho implies that K(fn) —)■ K{f) 
strongly in Ho- Hence, {fn, K{.fn))Ho (/> A'(/))-^p. Since A ■. H ^ H' is bounded, the 
mapping G :?{—)■ M, / i-A (A(/),/), is continuous. The linearity of A and property (iii) 
imply that G is also convex. Thus, G is weakly lower semicontinuous [5l Corollary 3.9]. 
Therefore, 

-W),f) = - {K{f),f)no < liminf ({A(/„),/„) - (A'(/„),/„)„,) = 0, 

n^oo 

because {L{fn),fn) —)■ 0 as n —?■ oo by assumption. We infer from hypothesis (i) that 
/ G -Af{L). But also / G A/'(L)-*-, so / = 0. Then, by hypothesis (iii), 

0 < Go = CoWfnfu < (A(/n),/n) = {K{fn),fn)Ho “ Wn),fn) ^ 0 , 
which is a contradiction. □ 

Let Ho = Lf. By la Prop. 2], assumption (ii) of Lemma ITO] holds. Furthermore, Lemma 
in] shows that (iii) holds true. Assumption (i) is a consequence of Lemma 0 Let f E V C H 
und set / = / — n^{f) G Af{L)-^. Then 

-(L(/), /) = (L(7), 7) > G||/||^ = G||/ - i77/)ll7 

since L{f) G Ll and {L{f),f) = {f,L{f))L 2 for f eV. This proves the hrst statement in 
Theorem El 


4. Explicit spectral gap estimate 

We present a second proof of the spectral-gap estimate dTOj) with explicit constants. The 
idea is to decompose the collision operator L into a mono-species and a multi-species part 
and to exploit the fact that the conservation properties of L are different from those of the 
mono-species part L™. Let assumptions (A1)-(A4) hold. 
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4.1. Decomposition. We decompose L = L™ + where L™- = (L™,...,L™), = 

and 

(31) ir(/<) = i«(/.,/.). i?(/)=X]^d(/../j)- 

Denoting by 77™' the orthogonal projection onto (with respect to the scalar product 

in Ll), we can decompose / according to 

(32) / = /ll+/a where/":= ;!”■(/), /!■:=/_/II. 

Lemma E] shows that 

(33) / G AA(L) if and only if // = + m • n + e|n|^) for Oj, e G M, u G 

(34) / G A/'(L™) if and only if // = M/^^(aj + Uj • n + ej|n|^) for Oj, Cj G M, M/ G 


and /II has clearly the form fl341) . 

For later use, we dehne the following bilinear forms 

(35) -(/, L^if))Li = 7 V / BuA,[h,]^M,M*dvdv*da, 

4 7r6xS2 

1 n „ 

(36) -(/, L '(/))/.2 = - V V / h,fM,M*dvdv*da, 

4 7r6xS2 

where hi = fi and 

Aj[7i] := h[ + 7'* - 7j - h*, Aij[hi, 7/ := 7' + h'* - hi- h*. 


4.2. Spectral-gap estimate for L™. Our starting point is the fact that the mono-species 
collision operator L™ has an explicitly computable spectral gap. A spectral-gap estimate 
for the linearized collision operator with n = 1 was proved in [26l Theorem 6.1, Remark 
1 ]: 

(37) ^ f BuAi[KfMiM:dvdv*da > ^ / (/, - n^{fi)fp,4v, 

4 7r6x§ 2 Poo,i 7 r3 

where = Am( 7 , Oi, O^) > 0, only depending on 7 , Oi, and (see (A3)-(A4)), can be 
computed explicitly, 

Di('y) := / Bii{\v — v*\, cos'd)M*dv*da, 

and i G {1,..., n} is hxed. This yields the following estimate for L™, where we recall that 
the space 77 is dehned in dH]). 

Lemma 11. With L™ defi'ned in (ITI]) . we have 

-(/,L™(/))i 2 > C^Wf - 77™(/)||^ for all f G Dom(L™), 
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where 
(38) 
and A™ 


cm ^ 




f3po 


A™'( 7 , Cl, C^) is given in (l37|) . 


Proof. We sum (ITtD over i = 1,... ,n and employ (l3^ to obtain 

n p 

(39) - (/, L”)/))!; > A” V / (/< - ^'“(/))"—<i«. 

i=l Poo,i 

It remains to estimate va in terms of z/j, defined in (jH]). The definition of Mj implies that 
Mj = {poa,j/Poo,i)Mi. This fact, as well as dehnition (|9]) of Oi, the lower bound ffTSj) . and 
assumption (A 6 ) give 


Vi = 




P. 


oo,j 


=lp° 


Bi.M*dv*da < V 
■ 


PooJ 


BiiM*dv*da = 

Poo,2 


We conclude that r'ii/poo,i > Vi/{{3poo), and inserting this bound into fl3^ yields the result. 

□ 


Lemma flTTll and the inequality — (/, L^(/))l« > 0 immediately show that 

-(/, L(/))^. > C^Wf - n^iml for all fev. 

However, we need the projection onto A/'(L)-*“ instead of A/'(L™’)-*-, which is contained in 
Af{L)-^. Therefore, we will exploit the part —{f,V’{f))Li to derive a sharper estimate. 

4.3. Absorption of the orthogonal parts. We prove that the contribution /-*- (intro¬ 
duced in (l32ld in the term —(/, L^(/))i 2 = + f^))Li can be absorbed by 

the Pi norm of /■*■. 

Lemma 12. Let rj = min{l, C^/S}, where C™' > 0 is given in Lemma 1711 Then, for all 

fev, 

-if, L(S))li > (c™ - 4,)||/ - f% - |(/l, L‘(/«))i;, 

where /II = i7™'(/) is the projection onto A/'(T™')-*-. 

Proof. By Lemma [TTl we hnd that 

(40) - (/, L(/))ij > C"‘|[/ - /« \\l - (/, L‘(/))i; > C'“||/ - /«||^ - r,(/, VU)),,. 

since —(1 — p){f, L'^{f))i 2 > 0 for 7 G (0,1]. We estimate hrst the expression Aij[hi, hf in 
dehnition (15B]) . writing hf = M~^^‘^ff and hf = 

Aij[hi, hjf = (Aij[hll, hj] + Aij[hf, hj]Y 

= + AAKi hff + 2Ay[hS, h}\\AiAhi, hf] 

> - AiAhf, hff. 
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Inserting this estimate into 


and 


gives 


n p 


BijAij[h}l, h''A^MiM*dvdv*da 


(41) 


*=i iA 


BijAij[h^, h-jYMiM*dvdv*da. 


We claim that the last term on the right-hand side can be estimated from below by 
IIup to a small factor. For this, we employ the invariance properties of Bij and the 
identity MiM* = M'-M'*'. 


BijAij[hY, hfYMiM*dvdv*da 


/R6xS2 


< 4 


< 16 


3x§2 


B^j{{{hiyY + {{hfy*f + {hyf + {{hfyf)M,M*dvdv*da 


Bij^kYY MiM* dvdv* da = 16 


5xS2 


/R6xS2 


Bij{fYYM*dvdv*da. 


Thns, the last term on the right-hand side of (HTD can be estimated as 


|EE 

1=1 jjti 


BijAij[hf, hfYMiM*dvdv*da 


'R6xS2 


n p Up 

> -477^^ / Bij{fYYM*dvdv*da >-47^^ / {fY^Uidv = -4:T]\\f 


±\\2 


taking into acconnt dehnition ([9]) of r'j. We infer from fl4TD that 


n p 

if,L(f))Li>{C"‘-mf-f'\\n + lY.Y. / 


Aij[hl, h YMiM*dvdv*da, 


and dehnition fl5B]) yields the conclnsion. 


□ 


4.4. Estimate for the remaining part. It remains to estimate the term —{fY L^(/ll)) 2 , 2 . 

Lemma 13. For /H G i.e. f} = mI^ {ai + Ui ■ V F eijnp) for some at, G M and 

Ui G we have 

f)b "• 
ti=i 

where > Q is defined in (l43ll . 


Proof. Thanks to the momentnm and energy conservation, we obtain differences of the 
momenta and energies, which will be crncial in the following: 

Ui ■ v' + Uj ■ v'* — Ui ■ V — Uj ■ V* = {Ui — Uj) ■ {v' — v), 
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|/| 2 | |/*|2 l |2 l *|2 / \ f \ l \2 l | 2 \ 

ei\v\ + ej\v \ — ei\v\ — ej\v \ = [ei — ej)[\v \ — |w| ). 

Using these identities in , h'!], where h- = ai + Ui ■ v + ei|np, we hnd that 




Ajj[hf, hh^MiM*dvdv*da 




=1 

B«(( 


M,- — M 


f) • {v' — n) + (cj — ej)(|U|^ — \v\^)YMiM*dvdv*da. 


Using the symmetry of Bij (thanks to assumption (Al)) and of MiM* with respect to v, 
the function G{v, v*, a) = Bij{ui — Uj) ■ {v' — n)(|n'p — |np) is odd with respect to {v, v*, a) 
and thus, the mixed term of the square in the above integral vanishes. Therefore, we obtain 

1 ”■ _ r 

-(/II,L'’(/Il))i 2 = Bij{\{ui-Uj) ■ {v' - v)\^ + {ei - ejY{\v'\'^ - \v\'^f) 

(42) 

Now, we claim that 


*=i jV 
X MiM*dvdv*da. 


'll . _ 'll . 

Bij{{ui—Uj)-{v' — v)yMiM*dvdv*da= * ^ 


Bij \v — v'\‘^ Mi M* dvdv* da. 


/M6xS2 " 3 Jr6xS2 

To prove this identity, we write Ui^k and Vk for the fcth component of the vectors Ui and 
V, respectively. The transformation {vk,vl,ak) i-A —{vk,vl,ak) for hxed k leaves Bij, Mi, 
and Mj unchanged but eA —such that 


/R6x§2 


Bijv'j^ViMiM*dvdv*da = 0 ioi i ^ k. 


Furthermore, 


/ BijVkViMiMj dvdv*da = 0 ioi i k, 

jRSxSS 

since the integrand is odd. Therefore, 


3x§2 


Bij{{ui — Uj) ■ {v' — v))"^MiM*dvdv*da 


^{ui^k-Uj^k){ui,£-Uj^e) / Bij{v'k-Vk){v[-vt)MiM*jdvdv*da 


k,l=l 
3 


/R6xS2 


= 'y\{ui,k - Uj^kf / Bij{vk - v',^)'^MiM*dvdv*da. 

k^i 7 r 6 x §2 

In fact, we can see that the integral is independent of k, and we infer that 


/R6xS2 


Bij{{ui — Uj) ■ {v' — v)YMiM*dvdv*da 
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k=l 




Bij\v — v'\‘^MiM*dvdv*da, 


from which the claim follows. 

Hence, can be estimated as 

3 

> T 5^ (|«i - (e, - e,)^). 

*J = 1 

where 

(43) = min [ Bij min I -\v — v'\'^,{\v'\‘^ — \v\'^y\ MiM*dvdv*da. 

Jjj6x§2 13 J 


It remains to show that > 0. The integrand of fl43l) vanishes if and only if |t'| = |t|. 
However, the set 

X = {(T,T*,cr) e X X : |u'| = |u|} 

is closed since it is the pre-image of {0} of the continuous function F{v, v*, a) = p, 

i.e. X = F“^({0}), recalling that v' depends on {v,v*,a) through ([2]). Since X ^ x 
X S^, its complement X^ is open and nonempty and thus has positive Lebesgue measure. 
Since the integrand in fl43|) is positive on X^, we infer that > 0. This finishes the 
proof. □ 


4.5. Estimate for the momentum and energy differences. The last step is to derive 
lower bounds for the differences ~ “ ^jY)- First, we recall some moment 

identities: 


(44) 



/ MiVjVkdv = pidjk, 
JR3 



Mi\v\‘^dv 


Ibpi 


for all 1 < i < n and 1 < j,k < 3. 


Lemma 14. Let f ^ with ff = + Ui- v -\- ej|f p) for 1 <i <n. Then 


fidv = pi{ai + 3ei), I fivdv = piUi, I fi\v\'^dv = pi{3ai + Ibci). 


A/2 


1 / 2 . 


Proof. Decomposing f = where /'I = n'^{f) and /"'■ = / — n"^{f), we infer from 

G (see fl34|) l that , ff-)L 2 = 0 and hence, by fl3T|) again, 

, f})L 2 = / Mi{ai + Ui-V + ei\v\Ydv = pi{ai + 3ei). 

7r3 

The other identities can be shown in a similar way. □ 


Lemma 15. For all f E V, we have 

n 

{\ui - Uj\^ + (ci - ejf) > 

*j=i 


- npml 
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where Ut, Ci are the coefficients of the ith component of in 

on Af{L), Ck > 0 is given by 


is the projection 


(45) 


Ck = 60 npoo max 

l<A:,£< 5 n 


/ 'i/JkfJinidv 


2 = 1 


and iffik) is an arbitrary orthonormal basis of J\f{L"^) in L^. 

Proof. We again decompose / = /'' + /"'■ with /H = i7™(/) and f^ = f — f^^- Then 

( 46 ) 11 / - ngmi < 2 {\\p\\i + ll/ll - 77 '-(/)iiy ^ 

We estimate first the difference g := — n^{f) = n"^{f) — n^{f) G ^{LC) (note that 

A/'(L) C Let {'ifk) be an arbitrary orthonormal basis of in Because of 

and VyVi G we have fjk G Pi. Then, by Young’s inequality, we find that 


,' 0 fc)L 2 V’fc 


k=l 


2=1 ^ 

5 n 

^ (^, 'fPk)Li {g, ifdbi (V’fc, i’eW 


5 n 


Ui(v)dv = '^{g,'lpk)Ll{,g,'ll^£)LlYl / '^kfJiVi{v)dv 


k,i=l 


2 = 1 


k,£=l 


5 n 


~ ’ “ k,i=l 

bn 


= 5n max \{'ipk,f^i)n\\\9\\h- 

Kk.Kbn ” Kk.Kbn 


k=l 


Thus, we infer from 


that 


11/ - n^U)\\u < 2||/^||^ + lOn max |(V^fc,||/ll - n^{f)\\l 2 . 

l<k,l<bn ” 

Because of Af{L) C Af{L'^), we have Tl'^n^ = 77^ and 

II/" - ^'■(/)Hi| = ll/"lli> - 2(/r”(/),77'-(/))i; + ||77y/)|||j 

= ll/'llij - 2(/,77y/))i. + ||77y/)||i5 = ||/l|||, - ||77y/)||'ij. 

Consequently, setting ko = lOn maxi<fc_£<„ |(//;fc, 

(47) 11/ - 77'-(/)||?, < 2||/-^||^ + A^„(||/"||i; - ||77y/)||i;). 

Next, we compute the norms of /II and n^{f). Moment identities flH|) show that 
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+ ISe^ + QaiCi). 


i=l 


For the computation of the norm of n^{f), we choose the following orthonormal basis 
(0j) = of Af{L) in Ll: 

0,, = = {^Poo)-^I‘^M\I\\v\^ - 3), 

where 1 < j < n and 1 < A; < 3. Then, using the moment identities of Lemma [T4l 


n+4 n 

I-^^(/)IIl 2 = = ^Pcx),i(“i + 3ej)^ + p, 

j=l i=l 


E 


2=1 
^M|2 


poo^l 

poo 


Hi 




Inserting the above identities for ||/''||^2 and ||i7^(/)||^2 into (l47jl . we conclude that 


11/ - ffp/)ll« < 211/ - -ff’”(/)||?, + ioPco E 


2=1 


Poo ,2 I I 2 
Poo 


Poo ,2 


E 7-“‘ 

2 = 1 


Poo,2 


6/coPoo ( ^ y-Ci 


i=i 


Then, if the inequalities 


(48) 

n 

\ '' Poo,i 1 |2 

E^i“-i - 

i=i 

n 

X Poo, 2 ^ 

(49) 

n / ^ 

\ ^ Poo,i 2 1 \ ^ Poo,/ 

/=1 Foo Y Poo 


*3=1 


*3=1 

hold, the lemma follows with Ck = 6/coPoo- 

It remains to prove fHHj) and fll9|) . To this end, we dehne the following scalar product 
on 


(u, v)p = ^ -^Ui - Vi, U= {ui, ...,Un),V = (ui, ...,Vn) G 
i=l 


nSn 


where Ui ■ Vi denotes the usual scalar product in The corresponding norm is ||m||p = 
{u,u)y‘^. Then 1 = (1,..., 1) G satishes ||l||p = 1. The elementary identity 

|2 n \2 _ ||„, ^ -1112 


Mp-{u,i)p = ||m- (m, i)pi|,p 


can be equivalently written as 

n 


2=1 


Poo,i I 12 

Poo 


Poo,i 


E 7-“‘ 


2 = 1 


E 

2 = 1 


Poo,2 
poo 


Ui 


-E 

i=i 


Poo p , 
Po 


-U, 


2 
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Then, using ELiPooj = Poo, 


j _ \ ^ Poo,i 


, Poo,I , 

1-I Ui 


E^(E 


2=1 

n 


p. 




Poo,2 
oo 




E 

, k^i 


Pc 


Poo,fc 
pot 


Poo,fc 
Pot 


E 


Poo J 

Poo 


u. 


E 

2=1 


Poo,2 
Poo 


-Mj) 




p. 


X]i3^i(Poo,j/Poo)('^i Uj) 


'Yhk^i Poo,A:/Pc 


^ ^ Aj (uj Mj ) 

j¥=i 


where = (Poo,i/Poo)(Efc^i(Poo,fc/Poo)) Since Y^j^i^j = 1, we may apply Jensen’s 
inequality to this convex combination, leading to 

2 


^^E^IE^l E^>i- 


Poo,A; 


Wo 


i=l 


P. 


Po 


, k^i 


Pc 






H Poo,i \ \ o PooJ I |2 ^ \ o I 

1-^ 1 2_^ - -\Ui -Uj\ < 2_^ \Ui - 


U 


Poo / • Poo 


'J I , 


*J=1 


since poo,j < Poo- This ends the proof. 
Now, we are able to prove Theorem [3l 


□ 


Proof of Theorem 0 By Lemmas [I2l [131 and UHl we obtain 

f,r)b EL 

-(/.n/))z,i > (C”-4r,)||/-/"||?( +ye (k - + (e. - eA) 

*J = 1 

The hrst term on the right-hand side is nonnegative if we choose rj = min{l, 4C'™'C'fc/(16C'fc-|- 
D^)}, and estimate ffTOll follows with A = rjD^/{8Ck)- □ 


5. Convergence to equilbrium 

In this section, we prove Theorem jH The idea of the proof is to adapt the hypocoercivity 
method of [29] to the multi-species setting. To this end, we need to verify the structural 
assumptions (H1)-(H3) in [29l Theorem 1.1]. The setting is as follows. 

Let L be a closed, densely dehned, and self-adjoint operator on Dom(L) C Ll such that 
L = K — A and the operators K and A satisfy the following assumptions: 

(HI) The operator A is coercive in the following sense: There exist a norm || ■ |Y on 
PL C Li and positive constants Pj (0 < i < 4) such that for all / G Dom(L) C "H, 

i>oll/lli|< >>.11/11?, <(/,A(/))i|<i> 2 ll/llt, 


(50) 
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(51) (V„/, V„A(/))ij > P3||V„/||?, - i>4f\\%. 

Moreover, there exists a constant Cl > 0 such that for all f, g ^ Doni(L), 

(52) {L{f),g)Li<CL\\f\\n\\9\\n- 

(H2) The operator K has a regularizing effect in the following sense: For all £ > 0, there 
exists C(£) > 0 such that for all / G Hi, 

{VJ,V,K{f))L 2 < e\\VJ\\l.+C{e)\\f\\l,. 

(H3) The operator L has a hnite-diniensional kernel and the following local spectral-gap 
assumption holds: There exists A > 0 such that for all / G Dom(L), 

-(/,^/))LS>A||/-i7^(/)||^, 

where H^ is the projection on 

Assumption (H3) is a consequence of Theorem [3l Next, we verify assumption (HI). 
Using Lemma [H] and the continuous embedding H ^ L^, we see that fISUD holds. For the 
proof of (1^ . we employ Young’s inequality: 


n « 

(V./,V.A(/))i 2 = V / v./.-v. 




E 

2=1 


> 


-E 

2 ^ 

2=1 


fi^vfi ' ^v^idv ~\~ 


Pi 


-fi dv + 


\Vyfi\‘^Uidv 

I'^vfil^cdv 


>^3||V./||^-P4||/||i., 

where z /3 = 1/2 and P 4 = maxi<j<„sup^gjjs |V^Pjp/(2Pj). Note that z /4 is hnite since 
is bounded and p* is strictly positive (see Lemma ED. Finally, inequality fl52D follows from 
the decomposition L = K — A, the compactness and hence continuity of K, the explicit 
expression for A, and the Cauchy-Schwarz inequality applied to {L{f),g)Li- 

It remains to verify assumption (H2). Let N := = (27r)“^/'‘exp(—| tP/ 4). We 

decompose K = it'd) _ where iFdi = {K'f\ ..., Kn^) and 


C = E 

i=i 


'M 3 x §2 


/' 


+ 


fi 


* M (M')V 2 (Mf)V 2 


dv*da, 


n « 

kP = V / B,,{M,M;Y/^f;dv*da 

^ iM3xS2 


for 1 < i < n. Because of M^M'^ = for all k, we hnd that 


K. 


“>(/) = E 

i=i ■ 


'R3x§2 




1/2 


M* 


(Mr)V 2 /' ^ (Mj)V 2 /^ 


(M;M'*)V 2 (MfM')V 2 


dv*da 
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i=i ' 

n 

E 

i=i ■ 




dv*da 


'R3xS2 




The transformation a i-A —a leaves v and v* unchanged and exchanges v' and v'*. Assump¬ 
tion (A5) {bij is an even function) ensures that Bij is unchanged under this transformation. 
Therefore, 


JR3x§2 

and we can write as 


Bijf'N'*N*dv*da = 


BijfyN'N*dv*da, 


(63) 

where 


if."' (/) = 6 E ( A (/<)+( A)). 


i=i 


KS'(Jk) = / S«(iV'*/; + N'jt)N’dv>da, 1 < s, j, i < 


n. 


/R3xS2 


Note that . In a similar way, we can decompose the operator 


K, 


( 2 ) 


n « n 

(/) = / ByfJN-dv'd^y = 5^(Poo,iP=oj)‘'''Jrf (/,), 

_ 1 J X j =1 


where 

K\f{fj) = N [ B,,N*f;dv*da. 

./r 3 x §2 

Next, we estimate the derivatives of Klj\ It is shown in [29l Eqs. (5.15)-(5.18)] that for 
all £ > 0, there exists C{e) >0 such that for any / G I < i, j, k < n, and £=1,2, 

(54) l|V„if'''(A)Hi, < e||V„A||i5 +C(e)||A||i,. 

Then we infer from (|53ll that 


n \ 




i=l j=l 


Ll 




*,j=l 


Ll 


< n(max ^ \\^vK^f {fi)\\l 2 . 
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Thus, by flSTj) . it follows that for £ = 1, 

n 

< n^(max Poo,i)^ (^ll + C(e)|l/i||i 2 ). 

V l<i<n ' ^ ^ ^ 

i=l 

A similar computation shows that this estimate also holds for 1 = 2. We infer that 

n 

||VA '(/)||^2 < 4n2(max Poo,i)^ 

V l<i<n ' ^ ^ V V/ 

i=l 

This proves assumption (H2) since e > 0 is arbitrary. 

Proof of Theorem\^ We have verihed that assumptions (H1)-(H3) are satished. Then, 
using exactly the same arguments as in the proof of Theorem 1.1 in |29], but now for the 
multi-species case, we conclude the exponential decay flT^ of the semigroup e*^, which is 
the hrst property of the theorem. 

It remains to show that the decay estimate ffTT|) follows from flT^ . For this, we write 
the initial value // as // = i7^(//) -l- (/ — i7^)(//), where is the projection onto Af{B) 
in Then the solution to (jS]) is given by 

fit) = e^^fj = e^^n^ifi) + e*^(J - n^)ifi), t > 0. 

We have already shown that 

||e‘B(/ - iT«)g||„._. < Ce-«||g||„._, for all g e Hl„, t > 0. 

In particular, the choice g = [I — i7^)(//) and the property (/ — 77^)^ = J — lead to 

||e*^(/ - n^)ifi)\\Hf^ < Ce-^^\\{I - n^)ifi)\\Hl^. 

It remains to prove that f^o = B^ifj) = e*^i7^(//) is the global equilibrium. Since 
BlJ^ifj) = 0 and U^ifj) does not depend on time, the constant-in-time function g = 
n^ifi) is the unique solution to the Cauchy problem 

dtg = Bg, t>0, git)) = n^ifj). 

This shows that n^{fj) = e^^g{0) = e^^U^ifj) and hnishes the proof. □ 
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